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Abstract 

The (planar) AdS Dirichlet problem has previously been shown to exhibit super- 
luminal hydrodynamic sound modes. This problem is defined by bulk gravitational 
dynamics with Dirichlet boundary conditions imposed on a rigid timelike cut-off sur- 
face. We undertake a careful examination of this set-up and argue that, in most cases, 
the propagation of information between points on the Dirichlet hypersurface is never- 
theless causal with respect to the induced light cones. In particular, the high-frequency 
dynamics is causal in this sense. There are however two exceptions and both involve 
boundary gravitons whose propagation is not constrained by the Einstein equations. 
These occur in i) AdS3, where the boundary gravitons generally do not respect the 
induced light cones on the boundary, and ii) Rindler space, where they are related 
to the infinite speed of sound in incompressible fluids. We discuss implications for 
the fluid/gravity correspondence with rigid Dirichlet boundaries and for the black hole 
membrane paradigm. 
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1 Introduction 

Because asymptotically AdS spacetimes are not globally hyperbolic, the associated gravita- 
tional dynamics requires boundary conditions at the timelike conformal boundary J? + (e.g. 
the Einstein Static Universe (ESU) for global AdS). At an operational level it is often con- 
venient to impose boundary conditions at a fixed timelike cut-off surface at some finite 
location and then to eventually remove the cut-off. However, it also of interest to consider in 
its own right the AdS Dirichlet problem defined by a finite boundary S^. For convenience we 
take Hd to be located at a fixed radial coordinate r = td, where r = oo is the AdS boundary. 
In principle one could keep either the region r > rr> (which we call the ultraviolet (UV) side) 
or the region r < rp (which we call the infrared (IR) side), though we will consider only the 
IR side below. In our contexts, the IR region will contain a smooth causal horizon as defined 
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by observers on E^. 

There are several motivations for this study ranging from the purely mathematical [1], 
to various conceptual issues involving AdS/CFT (see e.g. [2]), an understanding of the 
Wilsonian holographic renormalization group [3, 4], and finding an AdS/CFT embedding 
of the black hole membrane paradigm as inspired by [5, 6]. The long- wavelength limit of 
this problem was analyzed in some detail in [7], while earlier works [5, 6, 8, 9] discussed 
the analogous situation for Rindler spacetime. The structure of the resulting differential 
operators was addressed in [l]. 1 One might also be interested in more general boundary 
conditions on T, D (see e.g. [1, 3, 4, 9, 10]), but we focus on the Dirichlet problem below. 

Our present interest concerns an unusual feature found in [5, 7, 11]. These works studied 
the long- wavelength hydrodynamic regime in the context of a flat boundary Ed, either by 
computing linear quasi-normal modes [5] or [7, 11] by following the construction of the 
fluid/gravity correspondence [12]. 2 They found that when E^ was sufficiently close to a 
Killing horizon (td < r *n{T) for some critical r* D set by the temperature T) the effective 
speed of sound vs along E^ exceeded the speed c of null geodesies on E^; i.e., it became 
superluminal. 3 In fact this result follows directly from thermodynamics as in this regime 
the pressure fluctuations exceed the energy density gradients; i.e., dP/dE > c 2 . It is also 
consistent with the analysis of [5, 6] who, inspired by the black hole membrane paradigm, 
relate the gravitational dynamics of Rindler space with a cut-off to that of an incompressible 
Navier-Stokes fluid living on the cut-off surface (for higher order generalisations see [8]). 
Here we recall that an incompressible fluid has an effectively infinite sound speed and is 
invariant under a Galilean symmetry group. 

The natural concern is that transmission of information along the boundary at speed 
Vs would violate bulk causality. To make this connection tight, recall that in the linear 
regime the relevant solutions are small perturbations of the (planar) (d + l)-dimensional 

1 There it was pointed out that, even in the context of Ricmannian signature, the PDEs which determine 
the bulk geometry via Einstein's equations are not Fredholm with Dirichlet boundary conditions. This is 
due to a linearization instability (visible directly in the Gauss-Codacci equations or the radial Hamiltonian 
constraint) which arises near backgrounds for which Yijj has vanishing extrinsic curvature. We thank Toby 
Wiseman for explaining the result to us. As with more familiar linearization instabilities, we see no reason 
why this should call into question the self-consistency of the full non-linear problem. 

2 The analysis of [7] also relates the fluid living on So to the boundary fluid dynamics; this result will 
not play much of a role in our discussion. 

3 One way to gain intuition for this result is to note that the speed of sound is relatively independent of ru 
when measured using some fixed metric (associated with some fixed conformal frame) on the AdS boundary. 
The change in speed on comes primarily from transforming from fixed boundary coordinates to proper 
distance and time on Since the surface T,jj is infinitely redshifted in the limit where it approaches a 
horizon, the local sound speed on £c in this limit diverges. 
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Schwarzschild-AdS solution 



ds 2 = --f{r)dt 2 + 



r 



£ 2 dr 2 r 2 




f(r) = 1 




(1.1) 



or perhaps of Rindler space, which is a limiting case. Here dx 2 d _ x denotes the line element 
on (d — l)-dimensional Euclidean space, Causal curves satisfy ds 2 < 0, which for 

future-directed curves yields dt > f~ x \dx\. But since f(r) increases monotonically, curves 
on the IR side of have dt > f^ijo^dxl. The fastest causal curve between two points on 
T, D is thus a null geodesic on S fl and any faster transmission of information would violate 
bulk causality. 4 

One thus wishes to ask whether the above superluminal speed hints at a serious problem 
for Dirichlet boundary conditions. In pondering the issue, it is natural to compare with the 
behavior of a tachyonic (m 2 < 0) Klein-Gordon scalar in Minkowski space for which the 
group velocity v g = ^ = ^ k2 k +m2 diverges at k 2 = —m 2 > 0. Here it is well-known that the 
dynamics is nevertheless causal in the sense that it transmits no information faster than the 
speed of light. 

However, we caution the reader that the tachyonic scalar is not good model for our 
AdS Dirichlet problem. For the tachyonic scalar the divergence in v g is associated with an 
instability which manifests itself both at k 2 = —m 2 (where it turns on) and also at arbitrarily 
small k > (where v g is imaginary). In contrast, the AdS Dirichlet problem has been shown 
to exhibit stable hydrodynamic behavior in the limit of small k [5, 7, 11]; i.e., v g is mostly 
real with a small negative imaginary part. In fact, we see no reason to expect instabilities 
any k. 5 While a full proof of stability is beyond the scope of this work, we will comment 
further on the issue in §5. In particular, our work below will show that the system is stable 
in the large k limit and also provide numerical evidence for stability at intermediate k. 

The following logical possibilities remain: 1) the system somehow succeeds in violating 
bulk causality, 2) the previous analyses are subject to some subtle inconsistency, 6 or 3) 
the superluminal speed is merely an artifact of the gradient expansion. In the latter case, 
while the gradient expansion may still be useful for computing e.g., correlators on T^o in the 
limit of low wave number and frequency, it would approximate the effects of large-but-finite 

4 This is a property of the Schwarzschild-AdS background and our choice of a fiat boundary. In contrast, 
in global coordinates, there is typically a faster shortcut through the bulk when traveling between two points 
at the same r = ro, though so long as the null convergence condition holds these shortcuts disappear in the 
limit where Xd is taken to infinity [13]. 

See [14] for a discussion of why causal propagation with superluminal dispersion need not imply insta- 
bility. We will return to this point in §5. 

6 This option might be suggested by the peculiar features found by [7] when the Dirichlet surface was 
taken to be in the near horizon region. For instance the map between the boundary data and the Dirichlet 
data seemed to violate the gradient expansion, though they were not able to definitively pin-point the cause. 
However, we find this region to be amenable to a very clean analysis. 
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sources rather less well than in more familiar cases. In particular, despite the existence of a 
long- wavelength sound mode with v$ S> c, no matter how large and slow are the transmitter 
and receiver, if separated by a distance L the delay between emission and reception is always 
at least L/c. 

This short note is an attempt discern which of the above possibilities are realized. Were 
there no Dirichlet wall, the well-known fact that the Einstein equations are symmetric- 
hyperbolic would immediately exclude option 1 above. Unless the effect were somehow 
due to a boundary graviton, it is difficult to see how this conclusion could be modified 
by the presence of a boundary. But we can probe both this possibilty and also option 2 
above by analyzing the dispersion of linear perturbations at high frequencies. We expect 
the associated ultraviolet group velocity vuv = v g(k —> oo) along to be closely related 
to the true maximum speed of information transmission in the same way that this limit 
gives geometric optics propagation along null geodesies for Klein-Gordon scalars with any 
(positive or negative) m 2 . In most cases we find that vjjv precisely matches the speed c of 
null geodesies in arguing against option 1. The fact that the same equations lead to 
an infrared sound speed v$ > c and v uv < c sugests that the previous analyses are in fact 
consistent (i.e., that option 2 is not correct) and supports the idea that the superluminal 
speed found in [5, 7, 11] is an artifact of the gradient expansion (option 3). See e.g., [14] for 
a brief overview of (and references to) condensed matter systems with similar features. 

The problem of interest boils down to examining the quasinormal spectrum for spacetimes 
with a horizon and a rigid Dirichlet cut-off surface. We will undertake this analysis both for 
the well-studied Schwarzschild-AdS geometry and the Rindler spacetime. The latter arises 
as a near-horizon geometry of the former and in this near-horizon limit one encounters a 
divergent sound propagation velocity at low frequencies. The Rindler case is sufficiently 
special (and sufficiently simple) that we choose to treat it separately rather than carefully 
extracting the results from the more general Schwarzschild-AdS analysis. In particular, as 
discussed in section 5 the Rindler setting admits certain boundary gravitons for general 
boundary dimension d while the Schwarzschild-AdS case does so only for d = 2. 

The outline of this paper is as follows. In §2 we outline the general strategy for studying 
linear gravitational fluctuations both the Schwarzschild-AdS and Rindler backgrounds, re- 
viewing some of the standard material discussed earlier in the literature. As is well known, 
the fluctuation equations can be decomposed based on the transformation properties of the 
graviton modes under spatial rotation symmetries. We will argue that the relevant channel 
for our studies is one where the graviton fluctuations transform as scalars under the spatial 
rotation group. We then undertake a detailed analysis of the scalar channel equation in §3 
(for Rindler space) and §4 (for Schwarzschild-AdS), arguing that the long- wavelength su- 
perluminal sound mode smoothly evolves with wavelength into a causal, linearly dispersing, 
high-frequency mode. The analysis involves a combination of analytical work both to extract 
the precise long-wavelength dispersion and the UV behavior of the dispersion (using WKB 
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techniques), as well as some numerical studies to understand the intermediate wavelength 
dispersions. We conclude in §5 with a discussion of the implications of our result for the 
AdS Dirichlet problem and for the membrane paradigm. Some technical results relevant for 
our analysis are collected in the appendices. 



2 Gravitational fluctuations and Dirichlet boundary conditions 

The long-wavelength solution to the AdS Dirichlet problem considered in [7] begins by con- 
sidering fluctuations around the planar Schwarzschild-AdSd+i black hole (1.1), subject to 
Dirichlet boundary conditions at a fixed cut-off surface taken to be located at r = rjj. 
The focus there is on fluctuations which respect the Dirichlet boundary condition whilst 
remaining long- wavelength vis a vis the hypersurface metric. We take the metric on to 
be a fixed Minkowski metric on so as to discuss propagation of signals with respect to 

a fixed background light-cone (thus c = 1 henceforth). With the cut-off surface taken all the 
way to the AdS boundary, r D — > oo, the fluctuation equations for gravitational fields were 
first solved for Schwarzschild-AdSs in the seminal papers [15, 16]. 

We will examine a similar set-up, except that we are going to focus on the linearized prob- 
lem and allow for fluctuations with arbitrary frequencies/momenta. In order to understand 
causality issues, it will also prove to be useful to examine the behaviour of fluctuations in 
the Rindler geometry which arises as the near-horizon limit of the planar Schwarzschild-AdS 
black hole. 



2.1 Fluctuations of Schwarzschild-AdSs geometry 

The gravitational fluctuations can be decomposed following [17] into decoupled sectors based 
on the transformation properties under a spatial rotation group. In the radial gauge h Tll = 
with 

5G MN (t,x\r) = e ^ t+i " x h^(r) 5*5% z = x 1 , (2.1) 

we use the spatial rotation symmetry 5*0(2) transverse to z to decompose the modes into 
modes into scalars, vectors and tensors. Further, defining dimensionless frequency and mo- 
mentum variables 
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one finds the gauge invariant master fields in the three sectors to be [17]:' 



Tensor channel: = h 



X 

y 

u 



Vector channel: Zy = / _ rr (q h tx + Xoh z 



(ttT r 

11 

Scalar channel: Z s = (4 to q h tz + 2 ro 2 h zz + 2 q 2 h tt 



+ (h xx + h yy ) [q 2 (2 - f(u)) - ro 2 ] ) (2.3) 
Here we have introduced the useful dimensionless variable 

=► f{u) = l-u 2 (2.4) 

so that the domain between the horizon r = r + and the boundary r — > oo is mapped to a 
finite interval u£ (0, 1) with it = 1 being the horizon. 

In terms of these master fields the equations of motion are: 

Zt(u) - —— Z' T (u) + £± Z T (u) = 0, (2.5) 

uj uf 2 

Z v (u) + , /(q2/ _ m2) Z v {u) + — ^ Z v (u) = 0, (2.6) 



3r 2 (l+^ 2 ) + q 2 (2n 2 -3?i 4 -3) 

^(Srr^ + q^-S)) ZsM 

3ro 4 + q 4 (3-4M 2 + u 4 ) + q 2 (4M 5 -4u 3 + 4u 2 ro 2 -6ro 2 ) „ , . , „. 

+ i 2 / 2 — 5u = 0. (2.7) 

u f 2 (3 ro 2 + q 2 (m 2 — 3)) 

We are interested in solving the above equations subject to SGmn = on E d and also 
having a regular future horizon, which implies ingoing boundary conditions on various fields. 
In particular we require 

SGmn ->■ (l-u)~ iro/2 , as u ->■ I. (2.8) 

These equations simplify dramatically in the naive limit of large to and q limit with 
s = ^ = 0(1). Keeping only the formally leading terms in ro, all three sectors give rise to 
an universal equation 

Z"(u) + m 2 l -^4J~ z ( u ) = °- ( 2 - 9 ) 
uf 2 



7 The expressions for planar Schwarzschild-AdSd+i are given in [18]. We only discuss d — 4 here as it 
exemplifies the generic case. 
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This is obvious for the tensor and vector channels; for the scalar channel it follows from the 
fact that 3 + s 4 (3 - Au 2 + w 4 ) + s 2 (Au 2 - 6) = (1 - s 2 f) (3 + s 2 (u 2 - 3)) thus simplifying 
the last term in the equation. The universal nature of this limit reflects the fact that for 
J"d = oo thermal correlation functions in the dual CFT agree at high frequency with those 
of the CFT vacuum, and thus are controlled by conformal invariance. Since all channels 
correspond to operators of dimension A = 4, they behave identically in this limit. 

Applying a simple WKB argument to (2.9) immediately gives s 2 f = 1 on which 
is just the condition for propagation along the light cones of £j> To begin, note that 
WKB solutions to (2.9) are oscillatory near the horizon (where / = 0). Quasinormal modes 
require ingoing boundary conditions, so the magnitude \Z\ of the solution remains essentially 
constant until the WKB approximation breaks down at some vwkb where s 2 f = 1. The 
boundary condition Z{r D ) = thus requires r WKB < r D . But the case r WKB < r D may 
also be excluded since for r > twkb the WKB solution is a linear combination of rapidly 
growing and decaying real exponentials. Matching to the oscillatory solution requires the 
two branches to be of the same order in tt> at r = twkb- As a result, the solution Z(r) 
is exponentially large in ro at r = r B > t^/kb and cannot satisfy our Dirichlet boundary 
condition. 

We thus tentatively conclude twkb = ro, and thus s 2 f{rE>) = 1 as desired. Most of 
the rest of this work is devoted to showing consistency with the Rindler limit, setting up a 
framework in which one could compute sub-leading corrections, making the analysis of (2.9) 
more rigorous, and especially to showing that including the effects of poles in the naively 
sub-leading terms in (2.7) and (2.6) does not significantly change the result. 

2.2 The master equations in Rindler: near horizon analysis 

Another spacetime of interest for us will be the Rindler geometry whose metric we take to 
be given by 

o 2 

ds 2 = -%: dr 2 + dg 2 + dx t dx l (2.10) 

in d + 1 dimensions, where i is an arbitrary constant with dimensions of length. The above 
can be shown to be the near horizon limit of the planar Schwarzschild-AdSd+i black hole 
(see Appendix B). We might specialize if necessary to d — 4 but for the moment we can keep 
d arbitrary. The gravitational fluctuation equations are easy to obtain and the master fields 
follow from the results of [17] (see Appendix A). 

We take the cutoff surface in our Rindler space-time to be g = £, noting that the 
scale invariance of Rindler space makes all choices of I are equivalent. There is thus no 
loss of generality in setting the cut-off to be given by the scale I which secretly is the AdS 
scale before taking the near-horizon limit. Below, we denote the frequencies and momenta 
conjugate to d T and d x as w/£ and q/£ respectively. 
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For the fluctuations, we continue to work in the radial gauge h eil = with 

6G MN (T,x i ,Q) = e- i ? T+i i'hp,{Q)6&6% z = x 1 , (2.11) 

we use the spatial rotation symmetry SO (2) transverse to z to decompose the modes into 
modes into scalars, vectors and tensors. Note that regularity (an ingoing boundary condition) 
on the future horizon requires 

5G MN g- lw , as g^O. (2.12) 



The simplest channel to study is the tensor channel, where for modes with polarization 
in the spatial directions transverse to the momentum, one encounters a minimally coupled 
Klein-Gordon equation: 



V 2 Z 



T 







ld_ 
g dg 



[gZ' T (g)) + 



W" 



I 2 



Z t {q) = 0. 



(2.13) 



The vector channel equation can also be derived quite simply. One finds (with Zy defied 
as before in (2.3)): 



Z' v {q) + 



q 2 q 2 



w 



2 ()2 



g (w 2 £ 2 — q 2 q 2 ) 



Z' v {q) + 



w 



£ 2 



Z v {q) = 0. 



(2.14) 



In the sound channel we have an interesting curiosity: owing to the fact that a(r) — 1 for 
Rindler when compared with the general form given in (A.l), the master field degenerates 
to simply: Z s (g) = h(g) = Y^i^i hais)- We immediately see that the dynamical equation of 
motion reduces to 



Z'^g) + - Z' s (g) + 

g 



w 



I 2 



z s {q) = o 



(2.15) 



so that Zs(g) precisely satisfies the massless minimally coupled scalar wave equation as in 
(2.13). 

Applying the WKB argument of §2.1 to the present equations gives — — > 1 at high 
frequencies. As before, this tentative result matches the light cones on £/> It will be more 
rigorously established in §3 . 



3 Rindler Quasinormal modes 

Having derived the basic fluctuation equations of interest, let us now turn to a detailed 
analysis of the gravitational quasinormal modes. We begin with the Rindler geometry which 
as one might anticipate is amenable to exact analytic treatment (at least for the tensor and 
scalar channels). 



3.1 Quasinormal modes at low momenta 

Consider the Rindler spacetime (2.10); it is well-known that there are no quasinormal modes 
if we keep the entire region of the geometry including a part containing the Minkowski 
[19]. We are however interested in the region of spacetime within q = I which allows non- 
trivial modes satisfying (2.12) and vanishing on We focus on the scalar channel modes 
satisfying (2.15) (or equivalently the tensor channel modes satisfying (2.13)), which can be 
solved explicitly in terms of Bessel functions: 

Z s {q) = A J_ JW (-i q |) + B F_ iw (-i q |) , (3.1) 

where we take q positive and we have chosen the (arbitrary) sign in front of q to simplify 
treatment of a branch cut below. The ingoing boundary condition on the horizon sets B — > 0. 
Imposing Z$(q — i) — then leads to a transcendental equation for w 

q) = , or equivalently J_j w (— iq) = 0. (3.2) 



Since the zeros of the Bessel function are gapped, it is clear that there is no long-wavelength 
mode, which satisfies the above equation. These features are clearly illustrated in Fig. 1. 




-4 -2 (I 2 4 -4 —2 2 4 

Re(w) Re(w) 

q = q = 5 

Fig. 1: The quasinormal frequencies of the Rindler geometry in the scalar and tensor channels for two 
different values of the spatial momenta, q = (on the left) and q = 5 (on the right). The 
absence of low lying hydrodynamic quasinormal modes is clear in the low momentum limit. 
Because the generic expression degenerates at q = 0, the left plot was in fact generated using 
q= 1(T 9 . 
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The behavior we find for Rindler modes in the scalar and tensor channels is consistent 
with the linear analysis of [5] and the subsequent nonlinear analysis of [6] , who found that the 
Rindler geometry with a cut-off is equivalent as a dynamical system to an incompressible 
Navier-Stokes fluid. Now an incompressible fluid does not have any propagating sound 
mode, which is equivalently phrased in terms of there being no linearly dispersing mode in 
the system. 8 From the hydrodynamic analysis of [16] we expect that the linearly dispersing 
sound modes to appear in the scalar channel. Their absence in Rindler is a reassuring 
confirmation of the expected result. 



-\a - 



-1.5 - 



-2.0 - 



-2.5 - 




-.1 -2 -10 I 2 .1 

Re(w) 

Fig. 2: The quasinormal frequencies of the Rindler geometry in the vector channel for q = illustrating 
the low frequency shear mode. 

On the other hand one does expect to see a shear mode in Rindler geometry, owing to 
the presence of dissipative terms in the Navier-Stokes dynamics. This can be confirmed by 
numerical analysis of the vector mode equation (2.14), see Fig. 2. 

3.2 UV characteristics of Rindler quasinormal modes 

Having established that the Rindler geometry (2.10) with a cut-off at q = I does not have 
any linearly dispersing modes for q,w C 1 i.e., in the long-wavelength limit, let us now 
ask what the behavior of the modes is at high frequency. The reason for our interest is 
the following: although there is no finite-speed propagating long-wavelength sound mode 

8 We save comments on the pure pressure mode, which is in some sense a sound mode with strictly infinite 
speed, for §5. 
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in the Rindler geometry, the full system at long wavelengths is described [6] by a Galilean 
invariant Naiver-Stokes dynamics. Such dynamics clearly has propagation outside the light- 
cone of the hypersurface metric at g — £, which we have judiciously engineered to the be 
the Minkowski metric on M^ -1 ' 1 . The question we would like now to answer is whether the 
superluminal propagation associated with the hydrodynamic modes is a consequence of the 
long-wavelength approximation or whether it continues to badger the UV modes as well. 

To answer this question, we need to know the solution to (3.2) in the regime of q, w ^ 1. 
While the zeros of Bessel functions are easy to find numerically, this asymptotic regime is a 
bit more involved owing to the fact that we are simultaneously taking both the argument and 
the Bessel parameter to be large. Fortunately for us these asymptotics have been throughly 
worked out by Olver [20] (see also [21]). By a suitable use of the WKB approximation (a.k.a 
Liouville-Green approximation!) the large order double asymptotics of Bessel functions can 
be expressed in terms of Airy functions. In particular, for J v {vz) one has (for any positive 
integer n) 



J v {yz) 



i 

U3 



2k 



>fc=0 



V 



1 



z/ 2(n+l) 



i/3 



vfc=0 



V 



2k 



1 



2(n+l) 



with 



2.3 

- C 2 

3 S 



log 



i + vT 



-vr 



(3.3) 



(3.4) 



In (3.3), the functions Ak((), -Bfc(C) are polynomials of order 6k in C -1 ^ 2 an d (1 — z 2 )^ 1 ^ 2 
that are nevertheless finite at z = 1 (equivalently, ( = 0). In particular, the leading terms 
Aq, Bq at large v are finite constants. In (3.4) we follow [21] in placing the branch cut along 
the negative real z axis. 

For the case at hand from (3.2) we have 



l/=-!W, 



v z = — i q 



q 

w 



(3.5) 



implying that in the regime of interest (large u, fixed () the zeros j V)S of the Bessel function 
J v {z) are given in terms of those of the Airy function by 



Cs 



_ 2 
V 3 CL a 



(3.6) 



where a s are the zeros of Az(x). Recall that such zeros are real and negative (a s G K_) and are 
labeled by s G Z + . We also mention that the largest (least negative) zero is a\ ~ —2.33811. 
To avoid the branch cut in (3.4) for real q we use the symmetry of quasi-normal modes under 
w — > — w* to take Re(w) > below. 
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To find the UV dispersion relation we need only solve 



log 



w 



+ 7 



q J w 



2 1 

3 w 



(3.7) 



where we used the fact that the zeros of of Airy function are real and negative to write 
a s = e t7T \a s \. Since the r.h.s. has an overall w -1 it seems that for w, q ^> 1 one might forget 
about the r.h.s. and conclude from the l.h.s. that w ~ q. This is confirmed by a more 
careful analysis including the various phases. We thus conclude that for large w 



w ->■ v uv q 



with 



vuv 



(3.8) 



Re(w) 

q 



0.005 



0.004 



0.003 



0.002 



0.001 - 



logq 



2.5 x 10 4 



5x 10 4 



7.5 x 10 4 



10= 



Fig. 3: The asymptotic behavior of Rindler scalar (and tensor) channel quasinormal modes showing 
approach to the causal dispersion (3.8). 



The asymptotic dispersion relation (3.8) has corrections that will contribute both to the 
imaginary part of the dispersion relation and to the group velocity v g = . We will 

not try to extract them in detail; the interested reader may do so using the aforementioned 
asymptotics. One can however quickly check the approach towards the asymptotic behavior 
for instance by numerically solving (3.7) as shown in Fig. 3; the high momentum behavior is 
well approximated by 

w»q+ 1.856q£ + 0.115 q"s (3.9) 

In particular, note that v g approaches 1 from above as q — > oo. Similar behavior occurs 
for Schwarzschild-AdSs quasinormal modes; see Fig. 4 and Fig. 5. The imaginary part of w 
vanishes at leading order as may be seen by either by studying (3.7) directly by comparing 
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(3.7) with (3.4) and noting [21] that the latter maps real positive ( maps to real z < 1. In 
fact, we believe Im(w) may decay faster than any power law at large q. This belief is based 
on a brief investigation of the WKB approximation, though we leave a full analysis for future 
work. 



4 The Dirichlet Schwarzschild-AdSs quasinormal modes 

Quasinormal modes of Schwarzschild-AdSs are a well-studied topic. The beginnings were 
the early analysis of scalar fields in [22] and the seminal works of [15, 16] who discovered the 
connections to hydrodynamics. It is well known that when we impose Dirichlet boundary 
conditions for the fields Zs(u), Zt(u), Zy{u) on the boundary of Schwarzschild-AdSs, i.e., at 
u = 0, (td — > oo), we obtain hydrodynamic modes at low frequency. In particular, one finds 
dispersion relations: 

tt> = -J=q-0(q 2 ) , m,q«l (4.1) 

in the scalar channel [16], which implies a mode dispersing linearly with the conformal speed 
of sound vs = 4|- We are interested below in the behaviour of this mode as we move our 
cut-off surface Ed into the interior of the spacetime. 



4.1 Low momentum Dirichlet quasinormal modes 

One can ask how the sound mode behaves upon imposition of the Dirichlet boundary con- 
dition at r = td (equivalently u = up). The answer can be found either by analytically 
solving the scalar channel equation [5] in the regime of interest (tt>, q 1) in parallel with 
[17] or by using the fluid/gravity correspondence to obtain the non-linear conserved stress 
tensor on So and reading off the sound dispersion relation from its conservation [7] . Either 
approach yields the simple result 



ro = ±V _^ =, v o = ^_^. (4 . 2) 

In converting the dispersion relation between ro and q to a sound velocity we have accounted 
for the fact that the induced metric on the surface at fixed r = is the standard Minkowski 
metric once we rescale the t and x l coordinates by the local red-shift factors. The relative 
scaling is simply a factor of v/JItd) which has made its appearance in the formula for Vg. 
We further note that the sound velocity exceeds that of the local light speed on at small 
r [7]: 

Vg > 1 => «d> t r D <23r+ (4.3) 

V2 
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The prediction (4.2) compares very well with numerical solutions. We present the numerically- 
calculated dispersion relation Re(fo) as a function of q in Fig. 4 and Fig. 5 for various choices 
of ud- 9 

We wish to ask whether the superluminal propagation is restricted solely to the low 
momentum sector or does it also affect the high frequency UV behaviour of the theory. 10 
The asymptotic behaviour of quasinormal modes is hard to see directly from numerical 
studies owing to slow convergence towards the limiting value. However, we can study the 
asymptotic dispersion using a trick similar to that employed in §3. With this in mind we 
turn to developing a useful WKB approximation. 

4.2 A sound WKB for Schwarzschild-AdS 5 

For the reminder of this section we focus exclusively on the scalar channel master equation 
in Schwarzschild-AdS 5 (2.7) since it this channel that exhibits superluminal speeds at long 
wavelength. 11 The UV behavior in the shear and tensor channels is in fact easier to analyze 
as the WKB argument of section §2.1 applies without subtleties. In particular, the formally 
sub-leading terms in the tensor equation (2.5) are smooth bounded functions. In fact for the 
tensor channel a detailed WKB analysis was presented in [25] who analyzed the asymptotic 
behavior of quasinormal modes for a Klein-Gordon scalar in Schwarzschild-AdSd+i. They 
found that the large q behavior is given byro = q + 0(q~^) which was also confirmed 
numerically in [26]. Examining the vector equation (2.6), one does find a pole in the formally 
sub-leading first-derivative term. But this pole lies precisely at the same place where an 
analysis of the leading terms already indicates that the WKB approximation breaks down 

9 The numerical computation is done using the Frobenius method developed originally in [22]. This 
method whilst being quite straightforward, seems to suffer from poor convergence properties for large q as 
noted for instance in [18]. For this reason we restrict to q < 5. Nevertheless, a peculiar feature of the results 
is that the quasinormal frequencies seem to approach their asymptotic values rather slowly as q increases. 
This is quite clear in Fig. 4 and Fig. 5. While it would be interesting to understand this feature, we will 
will not study it in detail here. We comment only that it seems unlikely to be a numerical artifact. The 
point here is that, since we start our series expansion at the horizon, the Frobenius series solution should 
get better as Ud — > 1. Yet the gap between the numerically calculated tn at e.g., q = 5 and the asymptotics 
derived in §4.2 clearly increases with ujj. 

ln In [7] the authors argued that sensible long-wavelength dynamics on E^j can be obtained by 'projecting 
out' the sound mode. This involves taking the non-relativistic scaling limit described in [23, 24], and reduces 
the dynamics to a Galilean invariant incompressible Navier-Stokes system. However, this does not address 
the question of causality nor contain information about high frequency UV behaviour of the theory. 

n In fact, it is easy to verify that the shear and tensor channel equations have no low frequency linearly 
dispersing modes. While this has shown to be the case for the hypersurface E r> located at the AdS boundary, 
we have checked that this continues to hold by numerical studies of the equations (2.5) and (2.6). This is of 
course consistent with the non-linear analysis of [7]. 
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Fig. 4: The lowest quasinormal mode in the scalar channel at low and intermediate frequency and 
momenta for various choices of Dirichlet surfaces: ud = (top), up = 0.25 (bottom). We 
show both the real and imaginary parts of the dispersion relation. For Re(ro) vs q we indicate 
the prediction (4.2) for the sound velocity at low frequency using a dashed line, while the dotted 
line is the predicted value for the high-frequency modes (4.15). Note that the q = sound 
mode propagates causally for the above values of ud- As indicated in the text convergence to 
the asymptotic UV behaviour is in fact quite slow and hence we restrict to momenta q < 5. 
Finally, we also note that plots do not incorporate the red-shift scaling factor. 



and so does not affect the argument; i.e., this pole also leads to vjjv = 1- 

To develop a WKB analysis, let us start by writing the scalar channel equation in a 

12 The behavior of scalar quasinormal modes with large damping was obtained using WKB in [27]. More 
recently, geometric optics techniques have been used in moving mirror set-ups to analyze singularities of 
boundary correlation functions in [28, 29]; here however region between the boundary and the Dirichlet 
surface is of primary focus. 
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Fig. 5: The lowest quasinormal mode in the scalar channel at low and intermediate frequency and 
momenta for various choices of Dirichlet surfaces: U D = ^ (top), where the sound mode 
starts to propagate on the light-cone of the Dirichlet surface and ud = 0.9 (bottom) when 
the sound mode is acausal. The rest of the conventions are as in Fig. 4. 

Schrodinger form. We define 

ro = v q , A(u) = u 2 - u 2 pole , u vo[e = \Fiu v = a/3 (1 - v 2 ) , (4.4) 
and note that equation (2.7) is of the standard Strum-Liouville type: 

Z" s {u) + P{u) Z' s (u) + Q{u) Z s {u) = . (4.5) 

We may thus employ a wavefunction renomalization to convert the equation to Schrodinger 
form. Setting 

d P(v) 

Z s {u)=Z s {u)C{u) , -\ogC(u) = ^ (4.6) 

du 2 
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we obtain the wave equation: 



Z"(u) = (q 2 F(u) + g(u))Z(u), (4.7) 



with 



_ N l-u 2 -v 2 „, s 3-15u 2 + M 4 -13w 6 + t; 2 (9w 4 + 30u 2 -3) , , 

up{u) Au 2 f 2 (u) A(u) 

We would like to write down a WKB solution to (4.7) in the limit q ^> 1 taking into 
account two facts: (i) the function J-(u) has a zero at u v = ±\/l — v 2 and (ii) the function 
Q{u) has an additional pole at u — u po [ e . For the moment, let us ignore the latter problem 
and see how Olver's method sketched in [30] can be put to work. 

We acknowledge the fact that q ^> 1 and attempt to solve the equation (4.7) by dropping 
Q{u). This is essentially what we would do in standard WKB, but for the fact that we have a 
potential zero (hence a turning point) of J-{u) at u = u v within the domain of interest. Given 
this the domain of interest is u G (ud,1), we will assume that u v G (ud,1) for simplicity 
(later we will show that u v is at the edge of the domain i.e., u v = Ud)- Because of the 
potential zero, rather than directly writing a WKB ansatz it is useful to rewrite the equation 
(4.7) in terms of the Airy equation: y"(x) = xy(x). 

With this in mind, define a new variable ( such that ((u v ) = and ('(u v ) > so that 
near u = u v , ( ~ u — u v and further rescale the wavefunction Zs(u): 

^ ^ Z s =(f\ k W (4.9) 



du \ ( J \du 
so that the re-defined field solves 

d 2 W 



dC 

with 



(q 2 C + *(C))>V, (4.10) 



*(0 = + UFF'-dF 2 ) ^^ + C — • (4.11) 
vs; 16C 2 V ' 16 J 73 7 y 1 

Thus apart from the correction piece contained in \l/ we have reduced the equation into a 

known form at leading order in the q ^> 1. We then learn that the wavefunction to leading 

order in the WKB limit q 3> 1 (with v ~ 0{1)) is simply: 

Z a (u)=(jj 4 Ai(q 2 / 3 C) +£!(«), (4.12) 

where we have already ensured that the boundary condition (2.8) is satisfied. Given this, 
the boundary condition at u = Ud demands that 

Ai (q 2 / 3 Cd) = 0, (at u = u D ) Cd = ^ , (4.13) 
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with a s as before denoting Airy zeros. We can then finally conclude 



( D oc u d - u v ~ v = J 1 - u 2 D . (4.14) 

This is the desired result: the UV dispersion relation on the cut-off surface is consistent with 
the causal structure since the UV speed of propagation is given to be 

Vuv = T nc~ — 2~ = fr~ — 2~ = 1 • ( 4 - 15 ) 

q y/l-u 2 D 

It remains to work out the effect of the subleading term in the potential ty(u). We 
relegate this discussion to Appendix C. The key new feature is the pole at u — u po [ e in 
A, which a-priori precludes direct application of the theorem quoted in [30]. Nevertheless, 
Appendix C explicitly demonstrates that corrections to WKB analysis are small in the limit 
of large tt>. It follows that the UV dispersion relation (4.15) persists independent of the 
cut-off hyper surf ace's location td (or equivalently ud)- 

This indeed should have been expected from earlier analyses since it appears that the 
pole at Up i e is an artifact of the master field (e.g., compare the analysis of [16] and [17]). 
The individual equations of motion for the gravitational fluctuations do not contain the pole 
at all; they only fact regular singularities only at u = 0, ±1 and at u = ±v3. When we take 
linear combinations of the individual fluctuation equations to obtain the master field, new 
singularities can in principle arise. Using (2.3) the origin of u po i e appearing in the function 
A(w) can be traced back to the singularity at u = ±v3. 

Moreover, explicit numerical solution of (2.7) seems to be oblivious to the presence of the 
pole; e.g., the Frobenius expansion about u — 1 (the horizon) has a radius of convergence 
that reaches to the boundary at u = despite there being a regular singular point in the 
domain of interest [17]. The result may also be expected from the fact that the scaling 
limit described in Appendix B relates the Rindler geometry to the near-horizon region of the 
planar Schwarzschild-AdSs black hole, and from the fact that the Rindler sound channel is 
easy to analyze. 

In summary it appears that quasinormal modes of Schwarzschild-AdSs, which at low 
frequency are hydrodynamic sound modes with superluminal speed, smoothly morph into 
causally propagating modes in the large frequency limit. We also can confirm that in the 
shear and tensor channels the UV dispersion is linear with vuv = 1 as in (4.15). Thus, as 
anticipated in the beginning, the the asymptotic dispersion is precisely luminal. 



5 Discussion 

We have discussed the gravitational dynamics of the Dirichlet problem in AdS and in Rindler 
space in an attempt to understand the curious superluminal propagation of long- wavelength 
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modes seen in these examples in earlier analyses. By moving away from the long-wavelength 
regime, we were able to show that the high frequency (linear) gravitational fluctuations 
described by master fields in the bulk spacetime propagate causally with respect to the 
Dirichlet hypersurface. While we focussed on the scalar channel, we argued that the result 
extends to the vector and tensor channels as well. This follows from the universal behavior 
(2.9) of the master field equations of motion in the WKB limit. We expect similar results to 
hold for other bulk modes; e.g., scalar field modes would behave similar to the tensor modes 
of the graviton, while probe Maxwell fields will follow vector graviton modes. 

There remain several important issues to address which we discuss below: 

Boundary gravitons: The master fields used above are "gauge-invariant" in the sense 
that they are invariant under the action of all linearized diffeomorphisms. But such trans- 
formations can sometimes be physical, in the sense that the describe directions in phase space 
that are non-degenerate with respect to the symplectic structure. This typically occurs when 
the diffeomorphisms map to itself in some non-trivial fashion. In certain setting these 
transformations can even generate propagating physical modes called "boundary gravitons." 
The most famous example of this kind occurs for AdS3 with ro = oo [31]. 

Such boundary gravitons are not constrained by the Einstein equations and so a priori 
might propagate acausally. This acausality would be benign in the sense that, so long as 
the boundary conditions are invariant under diffeomorphisms, this symmetry will prevent 
the boundary gravitons from interacting with other degrees of freedom. But by considering 
more general (but still local) boundary conditions one would be able to construct systems 
in which essentially any source transmits information outside the light cones defined by the 
metric. 

It turns out that truly superluminal boundary gravitons appear in the AdS3 Dirichlet 
problem outside a BTZ black hole. Taking g to be proper distance in the radial direction we 
may write the metric in the form 

ds 2 = -f(g)dt 2 + h(g)dx 2 + dg 2 } (5.1) 

where we have considered the non-spinning case J = 0. 13 We normalize t by choosing 
J(Qd) = h(go) = 1 so that the induced metric on the hypersurface is the flat metric on 
M. 1 ' 1 . Consider fluctuations of the surface of the form g = gr> + e5g(t, x), with e< 1. To 
leading order the new surface has induced metric 

induced = -(1 + e f'(g D ) Sg(t, x)) dt 2 + (1 + e h\g D ) 5g(t, x)) dx 2 . (5.2) 

13 In the planar context one may obtain results for J ^ by acting with an appropriate boost in the 
x-direction. 
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The Dirichlet boundary conditions demands that (5.2) be flat. This yields 

{h'(g D )^ 2 -f'(g D )^j6g(t,x)^0, (5.3) 
implying that the boundary graviton mode 5r(t, x) propagates at all frequencies with a speed 



nap) 

v Dd=\hj7-w- (5.4) 



h'{Q D ) 



Since BTZ black hole geometries have, 



/ P Q _ r 2 p-Q \ 2 / p Q _i_ r 2 p-g \ 2 

f(g)=(- -i^) , Kq)=( +T f (5.5) 

g£?D _|_ ^.2 p — QD 

we find = i — > oo as go -> logr + . (5.6) 

Thus the boundary (Dirichlet) gravitons in BTZ are superluminal for all finite cut-offs and 
have a velocity that diverges relative to the speed of light on as this surface approaches 
the horizon. Since the dispersion relation is precisely linear, vdq is clearly the speed of 
information transmission. Hence this theory is truly acausal with respect to the light cones 
of the bulk metric. 

In contrast, the higher dimensional Schwarzschild-AdS geometries do not admit a bound- 
ary graviton. This may be seen by noting that, unlike the case of AdS3, the constraints 
imposed by the Gauss-Codacci equations suffice to kill all putative large diffeomorphisms. 
So for higher dimensional Schwarzschild-AdS the full dynamics is captured by our master 
field analysis and respects bulk causality. 

Curiously, however, the Rindler geometry admits an unusual boundary graviton in all 
spacetime dimensions. In particular, there is an extra mode with arbitrary time dependence 
so long as q is precisely zero. This mode exists because Rindler is a somewhat degenerate case 
for our Dirichlet boundary conditions. Recall that we fix the induced metric on the cut-off 
surface to be the flat Minkowski metric. Since the Rindler geometry (2.10) is itself a section 
of a higher dimensional flat space, we can ask on what co-dimension one surfaces of IR 1 '^ can 
we induce a flat metric. It is easy to see that this is satisfied on any surface g = g(r) in (2.10). 
Moreover, for such surfaces the extrinsic curvature is easily shown to have only K TT ^ 0. This 
then implies that radial canonical momentum ttmn also only has a vanishing rr component; 
ttmn oc 5 l M 5 3 N ] in short the radial momentum is a pure 'pressure mode' and does not change 
the boundary energy density as measured by the time-time component of the Brown- York 
stress tensor. This boundary graviton in Rindler has all the right properties to be the Rindler 
limit of the Schwarzschild-AdS sound mode. Recall that in the near-horizon Rindler limit 
we expect an incompressible fluid [5, 8, 9] and hence a sound velocity vs — > oo. This is 
consistent with the above behavior. In fact, one can identify the boundary graviton with 
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an enhanced symmetry of the Navier-Stokes equations in translationally invariant spatial 
domains. It corresponds to the freedom in non-relativistic Navier-Stokes equations to shift 
the pressure by an arbitrary function of time: P — > P + F(t). This symmetry is actually part 
of the enhanced symmetry of the Navier-Stokes equations, the Galilean Conformal Algebra 
[32]. It would be interesting to ascertain whether Dirichlet boundary conditions for Rindler 
actually lead to conserved charges associated with the full Galilean Conformal Algebra, and 
in particular whether the latter is the asymptotic symmetry group for the Dirichlet problem. 

As a result of this pure-pressure mode, the strict the Rindler Dirichlet problem is "acausal." 
But the same feature makes the dynamics ill-defined. Without an additional boundary condi- 
tion at spatial infinity, there is nothing to specify the time-dependence of this boundary gravi- 
ton. Regularizing the system by imposing a boundary condition on g(r) at large x restores 
causality, and we saw above that the same is true if we deform Rindler to Schwarzschild-AdS. 
We would expect similar features to hold for other physically-motivated deformations such 
as those discussed in [9]. 

Initial data and wave propagation: Modulo the somewhat bizarre behavior of AdS3 
geometries, we infer from the causal propagation of the UV modes that the bulk theory is 
indeed sensible with Dirichlet conditions imposed on a finite All that remains is for 
us to understand the implications of the acausal sound mode at low frequencies. This can 
be understood in simple terms as for example discussed in [14] in the context of optical 
filters. Consider a medium through which we have dispersive wave propagation with fre- 
quency dependent response. If the medium response is such that modes of lower frequency 
are damped relative to the higher frequency ones, one is automatically led to apparent super- 
luminal propagation. As is usually the case, causal propagation of modes can be translated 
into a statement about analytic properties of the response function. The analyticity con- 
straint is however sufficiently weak to allow superluminal group velocities even if there is no 
instability (tachyons, or amplification of modes). All that is really needed is a pattern of 
attenuation which is sufficiently asymmetric around the frequency of interest. To get such 
an effect in the long wavelength regime we require this to happen around zero frequency, 
which can be achieved for systems in which time-reversal symmetry is broken (e.g. by the 
sort of ingoing boundary conditions imposed in each example above). 

Stability: We have presented evidence that our systems are stable. At small momenta 
stability is implied by the results of [5, 7, 11], while at strictly infinite momentum it follows 
from our WKB analysis. The numerical results presented in Fig. 4 and Fig. 5 indicate stability 
at intermediate momenta. Furthermore, for the cases that reduce to the massless scalar 
wave equation, linear stability at all momenta follows from positivity of the energy and the 
existence of an everywhere-timelike Killing field outside the horizon. We remind the reader 
that sound channel in Rindler space is of this type. We leave a more complete analysis for 
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future work, though we note that the results of [33, 34, 18] may prove useful. 

Implications for the membrane paradigm: In cases with vs > c, it appears that 
the gradient expansion will correctly compute correlators on Ed in the limit of low wave 
number and frequency. Nevertheless, the expectation that bulk gravitational waves should 
propagate causally and the idea that the UV propagation speed v uv should better model 
the actual maximum speed of information transfer (as occurs e.g. for linear Klein-Gordon 
tachyons and in the same way that that the WKB approximation leads to the geometric 
optics approximation along null geodesies in a general spacetime) suggests that the gradient 
expansion approximates the effect of large-but-finite sources less well than in more familiar 
hydrodynamic settings. For example, sound waves produced by a human throat carry in- 
formation that propagates at the long-wavelength sound speed vs- But our results suggest 
that a source on our of any large but finite size transmits information at any arbitrarily 
slow bit-rate, this information nevertheless travels only at speed c < vg. 

Interestingly, whether vs > c or vs < c, the numerical results of figures Fig. 4 and Fig. 5 
suggest that the group velocity v g decreases with wavelength (and thus increases with q) 
at long wavelength. This suggests that long wavelength pulses tend to smooth out as they 
propagate, becoming less sharp. However, we have seen that when v$ > c, at sufficiently 
short wavelength the behavior changes and v g instead increases with wavelength so that in 
fact v uv < vir- Thus a pulse of compact support (in space or time) must in fact sharpen at 
the corners, and presumably flatten in the middle, as it propagates and damps. 

Finally, note that from the analysis in Appendix C we expect that tt> = q + (D(q~^. We 
expect the leading correction to have a (damping) imaginary part and thus to control the 
approach to equilibrium for high momentum modes. We note that this imaginary part is 
the same order in q as that found at high momentum in [25] using the scalar wave equation 
for the case r D = oo, suggesting that this UV equlibration rate is largely insensitive to the 
cut-off. The Rindler limit is special however. In this case the high frequency modes may 
take substantially longer to equilibrate, though we leave a detailed analysis for future work. 
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A The master fields and fluctuation equations 



Here we record the results of [17] concerning linearized gravitational fluctuations about 
backgrounds of the form: 



ds 2 = a{r) -f(r) dt 2 + dx * dx ' + & ( r ) ^ 



(a.i; 



i=l 



We reserve upper case Latin indices for the spacetime indices {£, x l , r} and lower-case Greek 
indices for indices tangential to a fixed r hypersurface {t, x % }. We impose the radial gauge 
h rfJ- = and consider fluctuations of the form (2.1). Using the remaining transverse spatial 
rotation symmetry SO(d — 2), [17] decomposes the modes into modes into scalars, vectors 
and tensors and find the gauge invariant fluctuations master fields to be 



Tensor channel: 
Vector channel: Zy 
Scalar channel: Zg 



a(r) 
1 



a(r) 
1 

a(r) 



(q h tx + xv h zx ) 
(2 to q hu + ttJ 2 h zz + q 2 h tt 



+ q 2 / 



1 + 



af 



a'f q 2 f 



h 



p — 1 



(A.2) 



B Schwarzschild-AdS — > Rindler: near horizon scaling regime 



Consider the Schwarzschild-AdSs metric given in (1.1). The proper distance to the horizon 
from a fixed radial position can be easily computed to be: 



dr _i 
V7~ 2 



cosh 1 



r = r, cosh - £ 



Armed with this information consider the coordinate change: 



r = r , cosh 



2s 



2rJ ' 



x = — X 



on the Schwarzschild-AdSs metric (1.1), which brings it to the form: 



ds 2 = e 2 



q tanh — g ) dr + cosh I — - ^> dx 



(B.l) 



(B.2) 



(B.3) 
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We are interested in this geometry with a cut-off imposed at 

g = t r 2 D = r\ cosh(2e) u = % = 1 . (B.4) 

r£, cosh(2£j 

While we will eventually focus on the limit e — > we note that the above metric is exact; 
we've just performed some diffeomorphisms on Schwarzschild-AdSs. 

Now let us take the limit e — > 0; in this limit, the geometry (B.3) simplifies significantly 
to become just 



ds = e 



dg 2 — %r dr 2 + dx 2 

I 2 



(B.5) 



which, up to an overall rescaling by a constant factor e 2 we recognize immediately to be the 
Rindler spacetime: 



2 



e 



which is the Rindler geometry. 

Before analyzing other aspects of the scaling introduced in (B.2), let us note that while 
(1.1) solves Einstein's equation with a cosmological constant A = — 1§ while the Rindler 
geometry (2.10) derived above, solves vacuum Einstein equations. This is perfectly sensible, 
since to get the Rindler spacetime we have to rescale by a scale factor e~ 2 as in (B.6). This 
scaling effectively rescales the AdSs cosmological scale to i = ie which vanishes to leading 
order in the scaling limit e — > 0. 

Note that the cut-off in Schwarzschild-AdSs at r = rp now implies the the Rindler 
spacetime comes equipped with a cut-off at g — I. As in the main text we denote the 
frequencies and momenta conjugate to the Rindler time translation and spatial translational 
symmetries, d T and <9 X as w/£ and q/£ respectively. 

To proceed, we also need a map between the frequencies and momenta frequencies and 
momenta in Rindler and those in the Schwarzschild-AdSs spacetime. To this end recall that 
in §2.1 we used dimensionless frequencies (2.2) 

w u(2 q ql 2 

m = 27r = 27;' q = ^f = 2^ + (R7) 

when we use T = Using the coordinate transformation and noting the extra factor of 
i in front of dr 2 we realize that the momenta in the metric (B.2) are related to those in 
Schwarzschild-AdSs via 

m = w, ^ = \\ (B.8) 

Consider now the near- horizon limit of the various equations in Schwarzschild-AdSs. 
Using (2.4) and (B.2) in the limit e — > 0, we learn 

dZ J 2 1 dZ 
du e 2 4 q dg 
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and that 

/(«)-> 4^. (B.10) 

This information then suffices to extract the leading e — > behavior of the Schwarzschild- 
AdSs equations of motion in the near horizon Rindler region. 

In the scaling region, the shear channel appears to be the worst (as one might expect 
owing to the non-trivial nature of the wave equation). The scalar and tensor channels are 
much simpler in the limit. The final expressions, after a bit of symbolic algebra is as given 
below for the three channels. 



Scalar/sound channel: The leading order terms in the e expansion (stripping off an 
overall 16fr -4) are: 



1 d 

g dg 



dZc 



dg 
20 g 2 



+ 



w 



Z s 



d 2 Z, 



+ 



3 dZ s 



3£ 2 I dg 2 5 g dg 



+ 



12 + 11 
\Qg 2 



7q 2 

10 1 2 



(b.ii; 



Tensor channel: The leading order terms in the e expansion (stripping off an overall 
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arc: 



1 d 

g dg 



dZ 



r 



+ e' 



dg 
20 g 2 



+ 



w 



t 2 



d 2 z 



7 9 dZx 
+ = r- + 



dg 2 5 g dg 



11 w 2 7q 2 

10 g 2 ~ 101 2 



(B.12) 



Shear channel: The leading order terms in the e expansion (stripping off an overall 



16§r pf) are: 



d 2 Z 



dg 



q 2 g 2 + w 2 £ 2 



+ £ 



2 20f? 2 



<i^ 2 



1 dZ 



v 



q 2 f) 2 ) q dg 



d 2 Zv 



+ 



w 



'V 



3q 4 ^ 4 - 16£ 2 q 2 w 2 ^ 2 +9w 4 f 1 rfZ 



ll£ 2 w 2 



10 g 2 



5(w 2 £ 2 -q 2 <^ 2 ) 2 
7q 2 
10 i 2 



g dg 



(B.13) 



We recognize that, in each equation, the leading term in e is indeed given by the Rindler 
fluctuation equations quoted in §2.2. There are further corrections at 0(e 2 ) (ignoring the 
overall factor we scaled out), which incoperate the corrections away form the near horizon 
region and contain both finite size corrections as well as the contributions from the (rescaled) 
cosmological constant . 
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C Corrections to the sound WKB for Schwarzschild-AdSs 



In this appendix we make good of our promise in §4.2 to estimate corrections to the WKB 
analysis. In particular, we want to be able to show that the corrections especially due to 
the pole are under control and that the conclusion regarding the UV dispersion relation 
continues to hold for asymptotic values of momenta. Following standard WKB tricks, let us 
write 

W(C) = Wo(C)e iaK) , (C.l) 
with = q 2 CVVo, i- e -> Wo(C) = Ai(q 2/3 (). Letting £ = q 2/3 ( correction is captured by 



dC 

a(£) satisfying 
with 



"(0 + 2 i a'(0 K(0 - (a'(£)) 2 = q~ 4/3 *(0 (C.2) 



m = W)~ e ' 2 - (c ' 3) 

da 



If ^(C) were bounded then we would have 1% ~ ty/K. However, we have an issue 

„2/3 ft 

from the pole of Q which we represent via ^(O = >_a p °' e locally. The pole dominates the 

S Spole 

correction a in the neighborhood of ( po i e , leading to ia = ( _T le ■ To see the cross-over let 

S spoZe 

us set up a perturbation theory for q ^ 1 and set 



«(0 = «o(0 + q ^ 7 i«i(0 + --- (C4) 

and write (C.2) as equations for ao and treat a\ as a perturbation. Consideration of the two 
limiting behaviors discussed above, suggests that we can truncate the equation for a to a 
linear system: 

iaZ(0 + 2ia f (0K(Z) = q-^*(0 . (C.5) 
We will post-facto justify dropping the (a' ) 2 term showing that it is (9(q _2//3 ) for consistency. 
Solving for «o i n (C.5) we formally find: 

ia (o = y , de/c- i (e) Jdiv&icd) , (c.e) 

with 

K(t)=exj>(2 [ k) ^elZ 3/2 , (C.7) 



where we naively estimate the integral of K(^). The remaining computation can be done 
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straightforwardly 

J J s Cpole 

Rpole f ,z -i(f 3 / 2 -f 3/ ? 



R 



,2/3 
pole 



(C. 



The latter integral can be estimated in terms of the exponential integral Ei(x): 



ia (0 



+ (1 + 4^) Ei (4^ e ( V / ^ - 
_ e 4 «pofe (2 + (1 + 4^ v^)) log(£ - W 



pole 



16q2/3£ 



pole 



4 + 2 lE (l + 2e 8 ^ (1 - 4e 



3/2 



•"pole 



*3/2 • 
pole- 



Ei 



3/2 - 
pole- 



+ (1 + 4(SM4U 



(C.9) 



which tells us that cto ~ (9(q~ 2 ^ 3 ) even including the contribution from the pole. In deriving 
the above we have only ascertained the effect of the pole and have ignored finite corrections 
from other contributions to the latter can all be shown to be (9(q -4 / 3 ). We also see 

the post-facto justification for dropping the non-linear term. All in all these results confirm 
our WKB analysis and indicate the the UV dispersion relation we derived is indeed robust. 
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